New insights into plant natriuretic peptide evolution: From the lysogenic conversion in Xanthomonas to the lateral transfer to the whitefly Bemisia tabaci.

Plant natriuretic peptide-like (PNP) are signaling molecules related to adaptive responses to stress. The Arabidopsis thaliana PNP (AtPNP-A) is capable of modulating catalase 2 (CAT2) and rubisco activase (RCA) activity in some circumstances. Interestingly, many plant-pathogens co-opted PNP-like molecules to their benefit. For instance, the citrus pathogen Xanthomonas citri carries a PNP-like (XacPNP) that can mimic and regulate plant homeostasis, and many phytopathogenic fungi carry effectors (e.g., Ave1 and AvrLm6) that are indeed PNP-like homologs. This work investigates the PNP-like evolution across the tree of life, revealing many parallel gains and duplications in plant and fungi kingdoms. All PNP-like proteins in the final dataset are structurally similar, containing the AtPNP-A active domains modulating CAT2 activity and RCA interaction. Comparative genomics evinced that XacPNP is a lysogenic conversion factor associated with a Myoviridae-like prophage identified in many Xanthomonas species. Surprisingly, a PNP-like homolog was identified in Bemisia tabaci, an important agricultural pest, being to date the second example of lateral gene transfer (LGT) from plant to the whitefly. Moreover, the Bemisia PNP-like homolog can also be considered a potential new effector of this phloem-feeding insect. Noteworthy, the whiteflies infest many plants carrying PNP-like copies and interact with some of their bacterial and fungal pathogens, strongly suggesting complex recipient/donor traits of PNP by LGT and bringing new insights into the evolution of host-pathogen arms race across the tree of life.

Introduction of SARS-CoV-2 C.37 (WHO VOI lambda) in the Sao Paulo State, Southeast Brazil.

The Lambda variants of interest (VOI) (C37/GR/452Q.V1/21G) was initially reported in Lima, Peru but has gained rapid dissemination through other Latin American countries. Nevertheless, the dissemination and molecular epidemiology of the Lambda VOI in Brazil is unknown apart from a single case report. In this respect, we characterized the circulation of the SARS-CoV-2 Lambda VOI (C37/GR/452Q.V1/21G) in Sao Paulo State, Brazil. From March to June 2021, we identified seven Lambda isolates in a set of approximately 8000 newly sequenced genomes of the Network for Pandemic Alert of Emerging SARS-CoV-2 variants from Sao Paulo State. Interestingly, in three of the positive patients, the Lambda VOI infection was probably related to a contact transmission. These individuals were fully vaccinated to COVID-19 and presented mild symptoms. The remaining positive for Lambda VOI individuals showed different levels of COVID-19 symptoms and one of them needed hospitalization (score 5, WHO). In our study, we present a low level of Lambda VOI circulation in the Sao Paulo State. This reinforces the essential role of molecular surveillance for the effective SARS-CoV-2 pandemic response, especially in regard to circulating variants.

Genomic monitoring unveil the early detection of the SARS-CoV-2 B.1.351 (beta) variant (20H/501Y.V2) in Brazil.

Sao Paulo State, currently experiences a second COVID-19 wave overwhelming the healthcare system. Due to the paucity of SARS-CoV-2 complete genome sequencing, we established a Network for Pandemic Alert of Emerging SARS-CoV-2 Variants to rapidly understand and monitor the spread of SARS-CoV-2 variants into the state. Through analysis of 210 SARS-CoV-2 complete genomes obtained from the largest regional health departments we identified cocirculation of multiple SARS-CoV-2 lineages such as B.1.1 (0.5%), B.1.1.28 (23.2%), B.1.1.7 (alpha variant, 6.2%), B.1.566 (1.4%), B.1.544 (0.5%), C.37 (0.5%) P.1 (gamma variant, 66.2%), and P.2 (zeta variant, 1.0%). Our analysis allowed also the detection, for the first time in Brazil, the South African B.1.351 (beta) variant of concern, B.1.351 (501Y.V2) (0.5%), characterized by the following mutations: ORF1ab: T265I, R724K, S1612L, K1655N, K3353R, SGF 3675_F3677del, P4715L, E5585D; spike: D80A, D215G, L242_L244del, A262D, K417N, E484K, N501Y, D614G, A701V, C1247F; ORF3a: Q57H, S171L, E: P71L; ORF7b: Y10F, N: T205I; ORF14: L52F. The most recent common ancestor of the identified strain was inferred to be mid-October to late December 2020. Our analysis demonstrated the P.1 lineage predominance and allowed the early detection of the South African strain for the first time in Brazil. We highlight the importance of SARS-CoV-2 active monitoring to ensure the rapid detection of potential variants for pandemic control and vaccination strategies. Highlights Identification of B.1.351 (beta) variant of concern in the Sao Paulo State. Dissemination of SARS-CoV-2 variants of concern and interest in the Sao Paulo State. Mutational Profile of the circulating variants of concern and interest.

A comparative genomic analysis of Xanthomonas arboricola pv. juglandis strains reveal hallmarks of mobile genetic elements in the adaptation and accelerated evolution of virulence.

Xanthomonas arboricola pv. juglandis (Xaj) is the most significant aboveground walnut bacterial pathogen. Disease management uses copper-based pesticides which induce pathogen resistance. We examined the genetic repertoire associated with adaptation and virulence evolution in Xaj. Comparative genomics of 32 Xaj strains reveal the possible acquisition and propagation of virulence factors via insertion sequences (IS). Fine-scale annotation revealed a Tn3 transposon (TnXaj417) encoding copper resistance genes acquired by horizontal gene transfer and associated with adaptation and tolerance to metal-based pesticides commonly used to manage pathogens in orchard ecosystems. Phylogenomic analysis reveals IS involvement in acquisition and diversification of type III effector proteins ranging from two to eight in non-pathogenic strains, 16 to 20 in pathogenic strains, besides six other putative effectors with a reduced identity degree found mostly among pathogenic strains. Yersiniabactin, xopK, xopAI, and antibiotic resistance genes are also located near ISs or inside genomic islands and structures resembling composite transposons.

Evidence of facultative parthenogenesis in three Neotropical pitviper species of the Bothrops atrox group.

We examined four suspected cases of facultative parthenogenesis in three species of a neotropical lineage of pitvipers of the Bothrops atrox group. Reproduction without mating was observed in captive females of B. atrox, B. moojeni and B. leucurus housed alone for seven years (the two former species) and nine years (the latter one). In addition to the observation of captivity data, we investigated molecularly this phenomenon using heterologous microsatellites. DNA was extracted from the mothers' scales or liver, from embryo and newborn fragments, and yolked ova. Four of the microsatellites showed good amplification using Polymerase Chain Reaction and informative band segregation patterns among each mother and respective offspring. Captivity information, litter characteristics (comparison of the number of newborns, embryos and yolked ova) and molecular data altogether agreed with facultative parthenogenesis predictions in at least three out of the four mothers studied: B. atrox (ID#933) was heterozygous for three out of the four markers, and the sons S1 and S2 were homozygous; B. moojeni (BUT86) was heterozygous for two out of four markers, offspring S1, S3, E2, and E4, and O1 to O6 were homozygous; and B. leucurus (MJJS503) was heterozygous for three out of four markers, and son E1 and O1 were homozygous. B. moojeni (BUT44) was homozygous for all loci analyzed in the mother and offspring, which although not informative is also consistent with parthenogenesis. This study represents the first molecular confirmation of different pitviper species undergoing facultative parthenogenesis among Neotropical endemic snakes.

Analyses of Seven New Genomes of Xanthomonas citri pv. aurantifolii Strains, Causative Agents of Citrus Canker B and C, Show a Reduced Repertoire of Pathogenicity-Related Genes.

Xanthomonas citri pv. aurantifolii pathotype B (XauB) and pathotype C (XauC) are the causative agents respectively of citrus canker B and C, diseases of citrus plants related to the better-known citrus canker A, caused by Xanthomonas citri pv. citri. The study of the genomes of strains of these related bacterial species has the potential to bring new understanding to the molecular basis of citrus canker as well as their evolutionary history. Up to now only one genome sequence of XauB and only one genome sequence of XauC have been available, both in draft status. Here we present two new genome sequences of XauB (both complete) and five new genome sequences of XauC (two complete). A phylogenomic analysis of these seven genome sequences along with 24 other related Xanthomonas genomes showed that there are two distinct and well-supported major clades, the XauB and XauC clade and the Xanthomonas citri pv. citri clade. An analysis of 62 Type III Secretion System effector genes showed that there are 42 effectors with variable presence/absence or pseudogene status among the 31 genomes analyzed. A comparative analysis of secretion-system and surface-structure genes showed that the XauB and XauC genomes lack several key genes in pathogenicity-related subsystems. These subsystems, the Types I and IV Secretion Systems, and the Type IV pilus, therefore emerge as important ones in helping explain the aggressiveness of the A type of citrus canker and the apparent dominance in the field of the corresponding strain over the B and C strains.

Population Structure and Genetic Diversity among Isolates of Coccidioides posadasii in Venezuela and Surrounding Regions.

Coccidioides posadasii is a pathogenic fungus that causes coccidioidomycosis in many arid regions of the Americas. One of these regions is bordered by the Caribbean Sea, and the surrounding landscape may play an important role in the dispersion of C. posadasii across South America through southeastern Mexico, Honduras, Guatemala, and Venezuela. Comparative phylogenomic analyses of C. posadasii reveal that clinical strains from Venezuela are genetically distinct from the North American populations found in (i) Arizona and (ii) Texas, Mexico, and the rest of South America (TX/MX/SA). We find evidence for admixture between the Venezuela and the North American populations of C. posadasii in Central America. Additionally, the proportion of Venezuelan alleles in the admixed population decreases as latitude (and distance from Venezuela) increases. Our results indicate that the population in Venezuela may have been subjected to a recent bottleneck and shows a strong population structure. This analysis provides insight into potential for Coccidioides spp. to invade new regions.IMPORTANCE Valley Fever is a fungal disease caused by two species of fungi: Coccidioides immitis and C. posadasii These fungi are found throughout the arid regions of North and South America; however, our understanding of genetic diversity and disease in South America is limited. In this report, we analyze 10 new genomes of Coccidioides posadasii from regions bordering the Caribbean Sea. We show that these populations are distinct and that isolates from Venezuela are likely a result of a recent bottleneck. These data point to patterns that might be observed when investigating recently established populations.

Origin and diversification of Xanthomonas citri subsp. citri pathotypes revealed by inclusive phylogenomic, dating, and biogeographic analyses.

BACKGROUND: Xanthomonas citri subsp. citri pathotypes cause bacterial citrus canker, being responsible for severe agricultural losses worldwide. The A pathotype has a broad host spectrum, while A* and Aw are more restricted both in hosts and in geography. Two previous phylogenomic studies led to contrasting well-supported clades for sequenced genomes of these pathotypes. No extensive biogeographical or divergence dating analytic approaches have been so far applied to available genomes. RESULTS: Based on a larger sampling of genomes than in previous studies (including six new genomes sequenced by our group, adding to a total of 95 genomes), phylogenomic analyses resulted in different resolutions, though overall indicating that A + AW is the most likely true clade. Our results suggest the high degree of recombination at some branches and the fast diversification of lineages are probable causes for this phylogenetic blurring effect. One of the genomes analyzed, X. campestris pv. durantae, was shown to be an A* strain; this strain has been reported to infect a plant of the family Verbenaceae, though there are no reports of any X. citri subsp. citri pathotypes infecting any plant outside the Citrus genus. Host reconstruction indicated the pathotype ancestor likely had plant hosts in the family Fabaceae, implying an ancient jump to the current Rutaceae hosts. Extensive dating analyses indicated that the origin of X. citri subsp. citri occurred more recently than the main phylogenetic splits of Citrus plants, suggesting dispersion rather than host-directed vicariance as the main driver of geographic expansion. An analysis of 120 pathogenic-related genes revealed pathotype-associated patterns of presence/absence. CONCLUSIONS: Our results provide novel insights into the evolutionary history of X. citri subsp. citri as well as a sound phylogenetic foundation for future evolutionary and genomic studies of its pathotypes.

Phylogenomics.

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. The abundance of genomic data for an enormous variety of organisms has enabled phylogenomic inference of many groups, and this has motivated the development of many computer programs implementing the associated methods. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Identification and analysis of seven effector protein families with different adaptive and evolutionary histories in plant-associated members of the Xanthomonadaceae.

The Xanthomonadaceae family consists of species of non-pathogenic and pathogenic γ-proteobacteria that infect different hosts, including humans and plants. In this study, we performed a comparative analysis using 69 fully sequenced genomes belonging to this family, with a focus on identifying proteins enriched in phytopathogens that could explain the lifestyle and the ability to infect plants. Using a computational approach, we identified seven phytopathogen-enriched protein families putatively secreted by type II secretory system: PheA (CM-sec), LipA/LesA, VirK, and four families involved in N-glycan degradation, NixE, NixF, NixL, and FucA1. In silico and phylogenetic analyses of these protein families revealed they all have orthologs in other phytopathogenic or symbiotic bacteria, and are involved in the modulation and evasion of the immune system. As a proof of concept, we performed a biochemical characterization of LipA from Xac306 and verified that the mutant strain lost most of its lipase and esterase activities and displayed reduced virulence in citrus. Since this study includes closely related organisms with distinct lifestyles and highlights proteins directly related to adaptation inside plant tissues, novel approaches might use these proteins as biotechnological targets for disease control, and contribute to our understanding of the coevolution of plant-associated bacteria.

High throughput techniques to reveal the molecular physiology and evolution of digestion in spiders.

BACKGROUND: Spiders are known for their predatory efficiency and for their high capacity of digesting relatively large prey. They do this by combining both extracorporeal and intracellular digestion. Whereas many high throughput ("-omics") techniques focus on biomolecules in spider venom, so far this approach has not yet been applied to investigate the protein composition of spider midgut diverticula (MD) and digestive fluid (DF). RESULTS: We here report on our investigations of both MD and DF of the spider Nephilingis (Nephilengys) cruentata through the use of next generation sequencing and shotgun proteomics. This shows that the DF is composed of a variety of hydrolases including peptidases, carbohydrases, lipases and nuclease, as well as of toxins and regulatory proteins. We detect 25 astacins in the DF. Phylogenetic analysis of the corresponding transcript(s) in Arachnida suggests that astacins have acquired an unprecedented role for extracorporeal digestion in Araneae, with different orthologs used by each family. The results of a comparative study of spiders in distinct physiological conditions allow us to propose some digestion mechanisms in this interesting animal taxon. CONCLUSION: All the high throughput data allowed the demonstration that DF is a secretion originating from the MD. We identified enzymes involved in the extracellular and intracellular phases of digestion. Besides that, data analyses show a large gene duplication event in Araneae digestive process evolution, mainly of astacin genes. We were also able to identify proteins expressed and translated in the digestive system, which until now had been exclusively associated to venom glands.

Worldwide Phylogenetic Distributions and Population Dynamics of the Genus Histoplasma.

BACKGROUND: Histoplasma capsulatum comprises a worldwide complex of saprobiotic fungi mainly found in nitrogen/phosphate (often bird guano) enriched soils. The microconidia of Histoplasma species may be inhaled by mammalian hosts, and is followed by a rapid conversion to yeast that can persist in host tissues causing histoplasmosis, a deep pulmonary/systemic mycosis. Histoplasma capsulatum sensu lato is a complex of at least eight clades geographically distributed as follows: Australia, Netherlands, Eurasia, North American classes 1 and 2 (NAm 1 and NAm 2), Latin American groups A and B (LAm A and LAm B) and Africa. With the exception of the Eurasian cluster, those clades are considered phylogenetic species. METHODOLOGY/PRINCIPAL FINDINGS: Increased Histoplasma sampling (n = 234) resulted in the revision of the phylogenetic distribution and population structure using 1,563 aligned nucleotides from four protein-coding regions. The LAm B clade appears to be divided into at least two highly supported clades, which are geographically restricted to either Colombia/Argentina or Brazil respectively. Moreover, a complex population genetic structure was identified within LAm A clade supporting multiple monophylogenetic species, which could be driven by rapid host or environmental adaptation (~0.5 MYA). We found two divergent clades, which include Latin American isolates (newly named as LAm A1 and LAm A2), harboring a cryptic cluster in association with bats. CONCLUSIONS/SIGNIFICANCE: At least six new phylogenetic species are proposed in the Histoplasma species complex supported by different phylogenetic and population genetics methods, comprising LAm A1, LAm A2, LAm B1, LAm B2, RJ and BAC-1 phylogenetic species. The genetic isolation of Histoplasma could be a result of differential dispersion potential of naturally infected bats and other mammals. In addition, the present study guides isolate selection for future population genomics and genome wide association studies in this important pathogen complex.

Morphogenetic characterisation, date of divergence, and evolutionary relationships of malaria vectors Anopheles cruzii and Anopheles homunculus.

The mosquito species Anopheles cruzii and Anopheles homunculus are co-occurring vectors for etiological agents of malaria in southeastern Brazil, a region known to be a major epidemic spot for malaria outside Amazon region. We sought to better understand the biology of these species in order to contribute to future control efforts by (1) improving species identification, which is complicated by the fact that the females are very similar, (2) investigating genetic composition and morphological differences between the species, (3) inferring their phylogenetic histories in comparison with those of other Anophelinae, and (4) dating the evolutionary divergence of the two species. To characterise the species we used wing geometry and mitochondrial cytochrome oxidase subunit I (COI) gene as morphological and genetic markers, respectively. We also used the genes white, 28S, ITS2, Cytb, and COI in our phylogenetic and dating analyses. A comparative analysis of wing thin-plate splines revealed species-specific wing venation patterns, and the species An. cruzii showed greater morphological diversity (8.74) than An. homunculus (5.58). Concerning the COI gene, An. cruzii was more polymorphic and also showed higher haplotype diversity than An. homunculus, with many rare haplotypes that were displayed by only a few specimens. Phylogenetic analyses revealed that all tree topologies converged and showed [Anopheles bellator+An. homunculus] and [Anopheles laneanus+An. cruzii] as sister clades. Diversification within the subgenus Kerteszia occurred 2-14.2millionyears ago. The landmark data associated with wing shape were consistent with the molecular phylogeny, indicating that this character can distinguish higher level phylogenetic relationships within the Anopheles group. Despite their morphological similarities and co-occurrence, An. cruzii and An. homunculus show consistent differences. Phylogenetic analysis revealed that the species are not sister-groups but species that recently diverged within the Kerteszia group, perhaps concomitantly with the radiation of bromeliads in South America or during the Pleistocene climate oscillations.

Biochemical, transcriptomic and proteomic analyses of digestion in the scorpion Tityus serrulatus: insights into function and evolution of digestion in an ancient arthropod.

Scorpions are among the oldest terrestrial arthropods and they have passed through small morphological changes during their evolutionary history on land. They are efficient predators capable of capturing and consuming large preys and due to envenomation these animals can become a human health challenge. Understanding the physiology of scorpions can not only lead to evolutionary insights but also is a crucial step in the development of control strategies. However, the digestive process in scorpions has been scarcely studied. In this work, we describe the combinatory use of next generation sequencing, proteomic analysis and biochemical assays in order to investigate the digestive process in the yellow scorpion Tityus serrulatus, mainly focusing in the initial protein digestion. The transcriptome generated database allowed the quantitative identification by mass spectrometry of different enzymes and proteins involved in digestion. All the results suggested that cysteine cathepsins play an important role in protein digestion. Two digestive cysteine cathepsins were isolated and characterized presenting acidic characteristics (pH optima and stability), zymogen conversion to the mature form after acidic activation and a cross-class inhibition by pepstatin. A more elucidative picture of the molecular mechanism of digestion in a scorpion was proposed based on our results from Tityus serrulatus. The midgut and midgut glands (MMG) are composed by secretory and digestive cells. In fasting animals, the secretory granules are ready for the next predation event, containing enzymes needed for alkaline extra-oral digestion which will compose the digestive fluid, such as trypsins, astacins and chitinase. The digestive vacuoles are filled with an acidic proteolytic cocktail to the intracellular digestion composed by cathepsins L, B, F, D and legumain. Other proteins as lipases, carbohydrases, ctenitoxins and a chitolectin with a perithrophin domain were also detected. Evolutionarily, a large gene duplication of cathepsin L occurred in Arachnida with the sequences from ticks being completely divergent from other arachnids probably due to the particular selective pressures over this group.

Biogeography and spatio-temporal diversification of Selenidera and Andigena toucans (Aves: Ramphastidae).

Andean uplift, Plio-Pleistocene climatic fluctuation, and river dynamics in the Amazon basin have all been implicated in the diversification of the South American avifauna. We reconstructed phylogenetic relationships in the genus Selenidera, which has served as a classic case of putative refugial speciation, and the closely related genus Andigena, to better understand the processes driving their diversification. Using mitochondrial and nuclear DNA sequences, we constructed a phylogeny to estimate the pattern and timing of divergence within and between seven lowland Selenidera toucanets and the five species of Andigena mountain-toucans, which together form a single clade. All phylogenetic analyses supported the monophyly of the montane genus Andigena, but indicated that the genus Selenidera is likely paraphyletic with respect to Andigena. Our time tree analysis is consistent with the orogenic uplift of the northern Andean range having initiated the divergence between Selenidera and Andigena, and that the movement and fragmentation of montane habitats in response to Pleistocene climatic oscillations likely influenced diversification within Andigena. Estimated divergence times for lowland Amazonian Selenidera did not support the Last Glacial Maximum (LGM) refuge hypothesis as an important biogeographic factor for the diversification of lineages studied here. The timing of divergence within Selenidera is consistent with the hypothesis that geographic isolation of areas of endemism generated by Amazonian river dynamics during the Plio-Pleistocene contributed to Selenidera speciation and current species distributions.

Synthesizing and databasing fossilcalibrations: divergencedating and beyond

Divergence dating studies, which combine temporal data from the fossil record with branch length data from molecular phylogenetic trees, represent a rapidly expanding approach to understanding the history of life. National Evolutionary Synthesis Center hosted the first Fossil Calibrations Working Group (3–6 March, 2011, Durham, NC, USA), bringing together palaeontologists, molecular evolutionists and bioinformatics experts to present perspectives from disciplines that generate, model and use fossil calibration data. Presentations and discussions focused on channels for interdisciplinary collaboration, best practices for justifying, reporting and using fossil calibrations and roadblocks to synthesis of palaeontological and molecular data. Bioinformatics solutions were proposed, with the primary objective being a new database for vetted fossil calibrations with linkages to existing resources, targeted for a 2012 launch.

Synthesizing and databasing fossil calibrations: divergence dating and beyond.

Divergence dating studies, which combine temporal data from the fossil record with branch length data from molecular phylogenetic trees, represent a rapidly expanding approach to understanding the history of life. National Evolutionary Synthesis Center hosted the first Fossil Calibrations Working Group (3-6 March, 2011, Durham, NC, USA), bringing together palaeontologists, molecular evolutionists and bioinformatics experts to present perspectives from disciplines that generate, model and use fossil calibration data. Presentations and discussions focused on channels for interdisciplinary collaboration, best practices for justifying, reporting and using fossil calibrations and roadblocks to synthesis of palaeontological and molecular data. Bioinformatics solutions were proposed, with the primary objective being a new database for vetted fossil calibrations with linkages to existing resources, targeted for a 2012 launch.

Temporal and spatial diversification of Pteroglossus araçaris (AVES: Ramphastidae) in the neotropics: constant rate of diversification does not support an increase in radiation during the Pleistocene.

We use the small-bodied toucan genus Pteroglossus to test hypotheses about diversification in the lowland Neotropics. We sequenced three mitochondrial genes and one nuclear intron from all Pteroglossus species and used these data to reconstruct phylogenetic trees based on maximum parsimony, maximum likelihood, and Bayesian analyses. These phylogenetic trees were used to make inferences regarding both the pattern and timing of diversification for the group. We used the uplift of the Talamanca highlands of Costa Rica and western Panama as a geologic calibration for estimating divergence times on the Pteroglossus tree and compared these results with a standard molecular clock calibration. Then, we used likelihood methods to model the rate of diversification. Based on our analyses, the onset of the Pteroglossus radiation predates the Pleistocene, which has been predicted to have played a pivotal role in diversification in the Amazon rainforest biota. We found a constant rate of diversification in Pteroglossus evolutionary history, and thus no support that events during the Pleistocene caused an increase in diversification. We compare our data to other avian phylogenies to better understand major biogeographic events in the Neotropics. These comparisons support recurring forest connections between the Amazonian and Atlantic forests, and the splitting of cis/trans Andean species after the final uplift of the Andes. At the subspecies level, there is evidence for reciprocal monophyly and groups are often separated by major rivers, demonstrating the important role of rivers in causing or maintaining divergence. Because some of the results presented here conflict with current taxonomy of Pteroglossus, new taxonomic arrangements are suggested.